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Abstract 

We consider an electron coupled to the quantized radiation field and 
subject to a slowly varying electrostatic potential. We establish that over 
sufficiently long times radiation effects are negligible and the dressed elec- 
tron is governed by an effective one-particle Hamiltonian. In the proof 
only a few generic properties of the full Pauli-Fierz Hamiltonian -ffpF en- 
ter. Most importantly, _Hpf must have an isolated ground state band 
for \p\ < p c < oo with p the total momentum and p c indicating that 
the ground state band may terminate. This structure demands a local 
approximation theorem, in the sense that the one-particle approximation 
holds until the semi-classical dynamics violates \p\ < p c . Within this 
framework we prove an abstract Hilbert space theorem which uses no ad- 
ditional information on the Hamiltonian away from the band of interest. 
Our result is applicable to other time-dependent semi-classical problems. 
We discuss semi-classical distributions for the effective one-particle dy- 
namics and show how they can be translated to the full dynamics by our 
results. 



1 Introduction 



Electrons, protons, and other elementary charged particles move in essence along 
classical orbits provided the external potentials have a slow variation. E.g. in 
accelerators it is safe to compute the orbits by means of classical mechanics. As 
a somewhat crude physical picture one imagines that the electron dressed with 
its photon cloud at any given instant t is in its state of lowest energy consistent 
with the current total momentum p(t). During the time span St the external 
forces change the momentum of the electron followed by a rapid adjustment 
of the photon cloud resulting in the new total momentum p{t + St). In this 
picture two physical mechanisms are interlaced. Since the forces are weak, the 
acceleration is small and radiative losses can be ignored. The effective dynamics 
of the dressed electron is conservative within a good approximation. In addition 
the total momentum p(t) has a slow variation and can be regarded as a semi- 
classical variable with respect to the effective Hamiltonian of the ground state 
band (lowest energy shell). The goal of our paper is to establish the validity of 
the physical picture. 

To set the scene, let us introduce the quantum Hamiltonian under consid- 
eration. In fact, as will be explained below, our technique is fairly general 
and uses only a few generic properties of the Hamiltonian. Still, it is useful 
to have a specific example in mind. We consider a free electron coupled to 
the quantized Maxwell field. The Hilbcrt space of states for the electron is 
H e i = L 2 (fii 3 ) and its time evolution is governed by the Hamiltonian — ^A. 
m is the mass of the electron and we have set % = 1. For the photon field we 
introduce the Fock space T\> = Hf over the one-particle space L 2 (IR 3 ) £g)(D 2 , 
i.e. T h = (&%L S n (L 2 (]R 3 ) <g)(D 2 )®" and S n the symmctrizer. Thus a state 
4> e T h is a sequence of vectors {<ft°\ ^\ . . .} with <j>W e 5„(i 2 (M 3 ) ®(E 2 )®" 
such that \\(f>\\ 2 = II 2 < 00 • On ^ b wc define the two component 

Bose field with annihilation operators a(k,X), where k £ IR 3 stands for the 
wave number and A = 1,2 for the helicity of the photon. The fields sat- 
isfy the canonical commutation relations [a(k,X),a*(k',X')} = S(k — k')S\\i, 
[a(k, X),a(k', A')] = = [a*(k, X),a*(k', A')]. The Hamiltonian of the free pho- 
ton field is the given by 



Hf=J2 I d 3 kLu(k)a*(k,X)a(k,X) 

\ i n 



(1) 



A=l,2 



with dispersion relation u)(k) = \k\, the velocity of light c = 1. The electron and 
the photon field are minimally coupled through the transverse vector potential 
A(x). To assure transversality we introduce the standard dreibein ei(fc), e2(k), 
k/\k\. Then 

A(x) - (2^)- 3 / 2 J d 3 k -J=e x (k) (e* k -*a(k,X) + e-* k -*a*(k,X)) . (2) 



A=l,2 ' 



A is an operator valued distribution. To make it an unbounded operator we 
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smoothen over the form factor p as 

A p {x) = [ d 3 x' p(x- x')A(x') (3) 



and assume that p is radial, smooth, of rapid decrease, and normalized as 
J d 3 x p{x) = 1. In the corresponding classical Hamiltonian p would be the 
rigid charge distribution. With all these preparations we can introduce the 
Pauli-Fierz operator of a free electron as 

H Q = — (-iV<g> 1 - eAJx)) 2 + 1 (g> H f (4) 
2m 

acting on 7i = 7i c \ ® H.t . e is the charge of the electron. In (0) x denotes the 
position operator of the electron on L 2 (H 3 ). Hq is invariant under translations 
jointly of the electron and the photons. Thus the total momentum 

P = Pel ® 1 + 1 <8> Pi , pi = ^ / d 3 k k a*(k, X)a(k, X) , (5) 

A=l,2 J 

is conserved, [Hq,p] = 0. This can be seen more directly by rewriting (Q) in 
momentum representation as 

1 2 
H = —(p e i®l-eA p (iV Pel )) +l®Hi (6) 

and then going to a representation in which p is diagonal. This is achieved 
through the unitary transformation T defined by 

n 

(T4>) {n) (p, fci, . . . , k n ) = 4> {n) { P (7) 

i=l 

where ip( n ' is the n-particle sector component of ip in electron momentum repre- 
sentation. The transformed Hamiltonian T~ 1 HqT 1 again denoted by Hq, takes 
the form 

H = ^(p-Pi~eA p (0)) 2 + H i . (8) 
We decompose Tt and Hq on the spectrum of p, 



H= d A pn p , H = d 6 pH (p). (9) 
Jtr 3 jtr 3 

The spaces H. p are isomorphic to Ti.{ in a natural sense and will be identified in 
the following. Hq(p) is just Hq acting on TLi for a given value of p. 

Physically one expects to have a dressed electron state for given momentum 
p, at least if \p\ < p c . In our semi-relativistic model states with \p\ > p c decay 
through Cherenkov radiation to lower momentum states. Thus provisionally we 
assume that there exists a p c such that for every p G A g = {p : \p\ < p c } H a (p) 
has a unique ground state with energy E(p), i.e. 

h (p)Mp) =e( p )Mp) (io) 
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has a solution ifjo{p) & Hf which is unique up to scalar multiples and E(p) = 
inf spec (iJo (??))■ We will discuss in Section 3 under which additional conditions 
on p our assumptions can be verified. 

E(p) is the ground state band energy. The corresponding projection is 
P g = J® d 3 pPo(p) where Pq(p) denotes the orthogonal projection onto the 
one-dimensional subspace spanned by ipo{p)- The states in RanP g are called 
dressed electron states. More explicitly 



RanP g = d^^Mp) G H; 4> £ i 2 (A g )| . (11) 
On RanPg we have 

V = ( d 3 p (e-< B WV(p)) MP) ■ (12) 



iH t„ 



Thus, if initially in RanP g , the dressed electron propagates like a single quantum 
particle with dispersion relation E(p) 1 which is generated through the interaction 
with the photons. In particular, it stays in the dressed electron subspace for all 
times. 

As already remarked, the motion of an electron is modified through external 
electro-magnetic potentials. In general, they have a slow variation on the scale 
set by the Compton wave length. Thus we add to Hq the external potential 
V(ex) (and possibly also an external vector potential A cx (ex)). e is a small 
dimensionless parameter which controls the variation of V. The external forces 
break the translation invariance of Hq and the total momentum is no longer 
conserved. A state ip initially in RanP g will no longer remain so under the time 
evolution generated by the full Hamiltonian P , which reflects that an accelerated 
charge looses energy through radiation. Since the external forces are weak of 
order £, we can expect radiation losses to be negligible. More precisely, the 
acceleration is of order e and by Larmor's formula the energy radiated over the 
time span r is £ 2 t. Thus the relevant time scale is of order e . On that time 
scale the radiation loss is of order e, whereas the cumulative effect of the forces 
is of order 1 and the electron moves on the scale set by the potential V. If the 
initial ip £ RanP g , the dressed electron should still be governed by an effective 
one-particle Hamiltonian, which is obtained from the dispersion relation E{p) 
through the Peierls substitution 

Pi =E{p-eA cx (ex)) + V(ex) : (13) 

and 

>- lHt ^ = I d 3 p (e- lHl V(P)) MP) + 0(e) (14) 



e 



for < t < e~ Y T with some suitable macroscopic time T. We will establish ( p"4| ) 
under some additional assumptions on Po and for A cx = 0. 

In (|l4|) it is crucial to assume that ip 6 RanP g . For a general ijj £ Ti. 
one expects that, in essence on a time scale of order 1, ip splits into outgoing 
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radiation and a piece which is approximately in RanP g . This second piece is 
then governed by the approximation (|l4|). To prove that this really happens is 
a problem of scattering theory for free electrons and outside the scope of our 
present investigation. To our knowledge this is an unsolved and very challenging 
problem. 

There are several difficulties with the picture proposed in (|l4|). Firstly the 
Pauli-Fierz Hamiltonian is infrared divergent. For p ^ 0, the photon cloud 
has an infinite number of photons, though of finite total energy. The physical 
ground state at fixed p ^ does not lie in Fock space. Even if we introduce a 
suitable infrared cut-off by assuming that p(k) — ► as k — ► 0, with p the Fourier 
transform of p, E(p) is not separated from the rest of the spectrum of Hq(p). 
At present, we have no technique available to control (|f4| ) in case Ho(p) has no 
spectral gap. To overcome both difficulties we are forced to give the photons a 
small mass, which means to set uj(k) — (m 2 h + fc 2 ) 1 / 2 . We emphasize that the 
radiation losses are not affected by this cut-off. 

As a second difficulty, which has been accounted for already, we observe that 
E(p) will cease to exist beyond a certain critical value p c , i.e. for \p\ > p c . This 
is most easily understood by considering the uncoupled Hamiltonian Hq(j>) at 
e = 0. ft has absolutely continuous spectrum and the only eigenvalue E(p) = 
2^-j-p 2 . This eigenvalue is isolated and below the continuum edge for \p\ < m e \ 
(m p h small) and is embedded in the continuum for \p\ > m c \. For small coupling 
the embedded eigenvalue should dissolve and E(p) exists only for \p\ < p c = m e h 

We could avoid the bounded extension of E(p) through a suitable modifi- 
cation of the boson dispersion u(k) and/or the electron dispersion p 2 /2m, cf. 
Section 2.3. However, the termination of bands is a fairly generic phenomenon. 
Thus, viewed in a more general setting, we have a classical Hamiltonian 

H cl (q,p) = E(p) + V(q) (15) 

corresponding to ( |l3|) on phase space T — H 3 x A g . The solution flow does 
not exist globally in time and for given initial conditions {q,p) there is a first 
time T (including T = oo) when the solution trajectory hits the boundary of 
T. This means we have to control the approximation (|1J) up to the time when 
substantial parts of the wave packet "leave" the allowed phase space. Beyond 
that time new physical phenomena appear not accounted for by ( [Ti]) . We believe 
that one of our main achievements is to develop a technical machinery which 
allows for such local approximations. 

It turns out that the validity of JlJ) relies only on a few rather general 
facts. Therefore we decided to prove (|l4|) as an abstract operator problem. The 
basic assumptions are the decomposition (^|) and a non-degenerate energy band 
of possibly finite extension separated by a gap from the rest of the spectrum. 
We mention three widely studied physical systems which posses this abstract 
mathematical structure, but are, at first sight, very different from the dressed 
electron. 

(i) For an electron moving in a periodic crystal potential, p becomes the quasi- 
momentum and the energy band is one particular Bloch band. Bloch bands are 
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discrete but may cross. The small parameter arises, as for Pauli-Fierz, through 
a, on the scale of the lattice spacing, slowly varying external potential. In |Q we 
studied the semi-classical motion in periodic potentials for Bloch bands that are 
isolated over the whole Brillouin zone. While we used a similar approach, the 
present paper is a substantial improvement in two respects. The approximation 
is local, no isolated bands have to be assumed, and we consistently avoid using 
any information on H(p) away from the band of interest. 

(ii) Electrons are lighter than nuclei by a factor of 2 • 10 3 -5 • 10 5 , which is the 
starting point of the Born-Oppenheimer approximation for molecular dynam- 
ics, p is now replaced by R, the collection of nucleonic coordinates. H(R) is the 
electron Hamiltonian for fixed nuclei and the band structure arises from eigen- 
values of H{R). They may cross or dive into the continuum. If we include the 
kinetic energy of the nuclei, the small parameter becomes e — {mci/m^c^s) 1 / 2 
and, except for the interchange of p and R, the Born-Oppenheimer approxima- 
tion fits our framework. Since in the following only bounded potentials will be 
covered — the kinetic energy is p 2 and plays the role of the potential — we post- 
pone a discussion of the time-dependent Born-Oppenheimer approximation to a 
separate paper. For an analysis from the point of view of wave packet dynamics 
we refer to || . 

(iii) The Dirac equation for a single particle has the electron and the positron 
band. One studies the motion of an electron, say, under slowly varying ex- 
ternal potentials [Q, |Q. As novel feature, the bands are doubly degenerate. 
Presumably this would also be the case if we include the electron spin in (JsJ) as 

H = — (a ■ (p - p f - eA p )f + H f . 
2m 

In the semi-classical limit the internal degrees of freedom (degeneracy) remain 
quantum mechanical and one has to approximate by matrix valued classical 
mechanics. 

Our paper is organized as follows. The abstract setting will be explained in 
Section 2. In this framework, almost by necessity, the theorems are stated as 
uniform convergence of certain unitary groups and of the corresponding time 
dependent semi-classical observables in the Heisenberg picture. The reader may 
worry, as we did, whether such convergence results imply the semi-classical 
approximation of quantities of physical interest, like position and momentum 
distributions. In fact, they do under very mild assumptions on the initial wave 
function. We could not find a coverage of sufficient generality in the litera- 
ture and therefore discuss semi-classical distributions in Section 6. Our main 
theorems are stated in Section 2 with proofs given in Sections 3 to 5. 

2 General setting and main results 

2.1 General setting 

As the "momentum space" M we take either lR d , d £ IN, in which case dp will 
denote Lebesgue measure on TR d , or a flat d-torus. In the latter case we take 
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dp to denote normalized Lebesgue measure on the torus. In the following Hf 
denotes any separable Hilbert space, although the notation should remind one 
of the Hilbert space for the Bose field in case of our main application. 

Let Hq be a self-adjoint operator on D(Hq) c H = L 2 (M) <g> Hf that can be 
decomposed on M as 

H = / d P H Q (p), (16) 

JM 

where {Ho(p),p £ M} is a family of self-adjoint operators with a common 
domain 2?o C Hf. We assume that the map p h- > Ho(p) is differentiable in the 
sense that for all p £ M and j = 1, . . . , d the limit 

(d P] H (p))(H (p) - i)^ := lim H ° {P ± = (g (p) - i)" 1 

exists in the norm of bounded operators on Hf. This defines, in particular, 
(V p Ho)(p) as a self-adjoint operator on Vq. 

For p in some compact and convex A g C M with non-empty interior let 
Ho(p) have an isolated non-degenerate eigenvalue E(p), i.e. 

h ( p )Mp) =e( p )Mp) 

with tpa(p) £ Hf and dist (E(p), a (H (p)) \ E(p)) > 0. We assume that E(-) £ 
C°°(Ag,lR). As before, Po(p) denotes the rank-one projection onto 4>o(p) and 
P g is defined as the orthogonal projection on RanP g given by (|ll|). The eigen- 
functions t^o(p), p £ A g , are defined only up to ap-dependent phase factor. One 



can choose this phase factor such that tpo(-) £ C (A g ,Hf) (cf. Lemma 11), and 
we assume such a choice in the following. On A g we require Hq(-) to be twice 
continuously differentiable in the same sense as above, i.e. for all p £ A g and 
j,k £ 1, . . . , d the limit 

{d Pk d Pj H (p))(H (p)-i) := hm — '- '- (H (p)-i) 

exists in the norm of bounded operators on Hf and depends continuously on p 
in the same topology. 

Let the potential V : JR d -> IR be such that V(x) = j dke lk ' x V(k), where 
J dk \k\ n \V{k)\ < oo for all n £ INo- This guarantees, in particular, that V and 
all its partial derivatives are in C°°(lR d ) and vanish at infinity. V e := V{ieS7 p ) 
is a bounded self-adjoint operator on L 2 (M) and the full Hamiltonian 

H = H Q + V E ® l n , (17) 

is self-adjoint on D(Hq). 

We define the corresponding one-particle Hamiltonian on L 2 (M) by 

Hf = E(p) + V E . (18) 

Note that E(p) is a priori only defined for p £ A g . To make things simple, we 
continue E(p) to a smooth and compactly supported but otherwise arbitrary 
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function on M. Since we will be interested only in the behavior for p e A g , 
we use this global function to define H\ . The corresponding unitary groups are 
denoted by U(t) = e~ im and U x {t) = e - lHlt . 
Our goal is to show that, in a suitable sense, 

(U(t/e)-U 1 (t/e))P g = 0(e) (19) 

for macroscopic times t < T < oo, where T depends only on V and the initial 
momenta. 

But before we can state the precise result, we have to make sense of Hi 
acting on RanP g . According to (|Tl]) we define the map U : RanP g — > L 2 (M) 



U(<Wo) = 4>. (20) 
Its adjoint U* : L 2 (M) — > RanP g is given by 

U*<\> = J S ' d P K As ( P )<f>(p)Mp), (21) 

where here and in the following 1a denotes the characteristic function of a set 
A. Clearly U is an isometry and U*U = 1 on RanP g . 

The effective dynamics cease to make sense, once the momentum of the par- 
ticle leaves A g . If that is not excluded by energy conservation, the comparison 
( pl)| ) can hold only up to some finite time, which can be determined from the 
classical dynamics generated by 

H cl (x,p) = E(p)+V(x) 

on phase space H d x M. Let A; C M, the "set of initial momenta" , and A m C M 
some "maximal set of allowed momenta" , both be compact. For AcM compact 
and 5 > we let 



A + 5 := pe M : inf \p - k\ < 8 
{ feeA 1 ~ 

be the corresponding <5-enlarged set, which is again compact. With : IR d x 
M — > M denoting the momentum component of the classical flow, we define for 
Ai + S c A m 

l£(Ai, A m ) := sup {t : supp ((l Ai o % s )(x, .))+iC A m Vs e [0,t],xG M d } 

as the maximal time for which the momentum of any (i.e. with any starting 
position) classical particle with initial momentum in A; stays within a 5-margin 
inside of A m . 

The following lemma shows that the classical bound on the momentum is 
respected also by the quantum dynamics in the limit e — > 0. Let P; = 1ai(p) 
and P m = t Am (p). 
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Lemma 1. Let Ai, A m C M be both compact and A; + 5 C A m for some 
5 > 0. For any T < oo with T < (A ; , A m ) there is a C < oo such that for all 

te[0,r] 

't' 



(1 - Pn) t/l - P 



(22) 



£(L2(M)) 



For the proof of Lemma [j] see Section 3. 
2.2 Main results 

In the following we will always assume that A; C A g and we can therefore 
identify P; with U*PiUP g to keep notation simple. I.e., Pi projects on "dressed 
electron states" with momenta in Ai. 



Theorem 2. Let A; be compact with Ai + i5 C A g for some <5 > 0. For any 
T < oo with T < T£(Ai, A g ) there is a C < oo such that for all t G [0, T] 



U[-} -U*U X - M Pi 



< Ce. 



(23) 



C(H) 



This means, in the case of our main application, that dressed electron states 
with initial momenta in Ai evolve according to the effective dynamics without 
radiation losses for times of order e~ l as long as the momenta of the correspond- 
ing classical orbits stay inside A g . Note, in particular, that the macroscopic time 

depends only on Ai, but not on e. 

The equivalence between the full dynamics and the effective one-particle 
dynamics on the isolated energy band extends to the level of semi-classical or 
macroscopic observables. On a suitable domain in 7i the macroscopic (Heisen- 
berg) position operator is given by 

x £ {t) = U (-t/e) {eiVp <g> 1) U {t/e) 

and its effective one-particle counterpart on H 1 (JR d ) by 

x\(t) = Ui {-t/e) eiV p Ui {t/e) . 



Theorem 3. Let Ai be compact with Ai + S C A g for some 5 > 0. For every 
T < oo with T < T£(Ai, A g ) there is a C < oo such that for all t € [0, T] 



(x e {t)-U*x\{t)U 



P 



C(H) 



< Ce. 



(24) 



Theorem [| is not a trivial corollary of Theorem g for several reasons. It is a 
statement about unbounded operators, whose difference, however, is bounded. 
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More importantly, the difference of the position operators does not vanish even 
at time t — because [a; e (0),P g ] ^ 0. 

More generally, consider classical symbols a S C°°(JRr x M) such that for 
all multi-indices a, (3 there exists C a .p < oo with 

sup \d£dPa(x,p)\ < C a ,p ■ 

x.p 

We use the notation of Q and denote this set of symbols by 5q(1). For a € <Sq(1) 
Weyl quantization leads to an operator a W ' £ S C(L 2 (M)) given by 

(a w ' e 0) (p) = (2^)~ d J dxdka (ex, e-^ h > x 0(A) , (25) 

with operator norm that is bounded uniformly in e (cf. Section 4 for details). 
For a € Sg(l) let 

a e (t) = [/ (_i/ e ) ( a w ' £ ® l) 17 (i/e) (26) 

and 

of (t) - f/i (-t/e) a w < £ Ui (t/s) . (27) 



Theorem 4. Let A; be compact with A; + S C A g for some S > and 
a € 5g(l) n C 0O (]R d x M). For every T < oo with T < r*(Ai, A g ) there is a 
C < oo such that for all t e [0, T] 



(a £ (t) -Wali^U^Pi 



C(H) 



< Ce. 



(28) 



Coo denotes the set of continuous functions vanishing at infinity. Presumably 
Theorem [I] holds as well for a G £o(l), however, the assumption a £ Coo allows 
for a simpler approximation argument in the proof. 

2.3 The massive Nelson and Pauli-Fierz model 

It is of interest to see whether our abstract assumptions are in fact satisfied for 
physically, at least semi-realistic models. The best studied case is the Nelson 
model H , where the coupling is to the position of the particle and the Maxwell 
field is replaced by a scalar field. Switching immediately to the total momentum 
representation, cf. (||), the Nelson Hamiltonian reads 

H^{p)= l -{ P - Pi f+H l + {2n)- d ' 2 J dk-^={a{k)+a*(-k)) , (29) 

where instead of (Q), (pi), a(k), a*(k) is a one-component Bose field over lR d . 
Again, p is the total momentum and regarded as a parameter. If J dk |p| 2 (w~ 3 + 
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u> 1 ) < cxd, then H^(p) is bounded from below and self-adjoint with domain 
D{H t ). 

According to a result of Frohlich under this condition with w{k) = (m 2 h + 
fc 2 ) 1 / 2 , m p h > and in dimension d = 3, i?N has an isolated, nondegenerate 
ground state band for \p\ < p c withp c > V%— 1. If in (|2^) we replace the electron 
dispersion by E (p) = (m^+p 2 ) 1 / 2 , m e i > 0, still keeping = (m^ + A; 2 ) 1 / 2 , 
Tiph > 0, then at zero coupling the ground state lies strictly below the continuum 
edge for all p, which means p c = oo for p = 0. As proved in pj p c — oo persists 
to arbitrary coupling strength. 

A larger class of bosonic dispersion relations is studied in jT3. For the 
particular case Lu(k) = luq > and dimensions d = 1, 2 one has p c — oo, whereas 
for d = 3 and small coupling indeed p c < oo J§| . 

Frohlich encloses the Bose field in a periodic box and proves that as the box 
size goes to infinity, there is still spectrum of (p) strictly below the continuum 
edge, provided \p\ < p c which has to be estimated through a separate argument. 
Thus Hn{p) has a ground state for \p\ < p c . To establish uniqueness he uses 
that e^*^^ in Fock space representation is positivity improving, except for a 
small error which goes to zero as t — > oo. From analytic perturbation theory 
one concludes then that the ground state band is isolated and analytic in p, 
provided \p\ < p c . Hence, the massive Nelson model satisfies the assumptions 
underlying Theorem |2[ [}] and 0. 

For the Pauli-Fierz model (H) the existence part of the proof goes through 
without essential changes. Unfortunately, the positivity part fails and no ar- 
gument is known to ensure the uniqueness of the ground state for |p| < p c . 
Since for e = and small p the ground state band is isolated, it must remain 
so for small e by perturbation theory, however, with the undesirable feature 
that smallness depends now on m p h. Even for the massive Pauli-Fierz model a 
general proof of uniqueness would be most welcome. 

3 Preliminaries from semi-classical analysis 

Although Theorem]^ compares a full quantum evolution to an effective quantum 
evolution and says nothing about semi-classics, its proof makes use of Lemma 
|l|. To prove Lemma |l| we need some tools from semi-classical analysis, which 
are introduced and applied to the one-particle Hamiltonian Hi in this section. 
Furthermore, in Section 6 we will consider semi-classical distributions based on 
the notions and results discussed here. 

Since E and V are both smooth bounded functions, we can apply standard 
results of semi-classics to the Hamiltonian 

flifofeVp) =E{p) + V(ieV p ) 

acting on L 2 (M), where the roles of momentum and position are exchanged and 
the role of h is taken by e. In the following we will simply ignore this difference 
and leave the necessary changes as compared to the standard case to the reader. 
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Note, however, that the change in sign, iV instead of — iV, is canceled by the 
fact that p' = q and g' = —p is the canonical transformation interchanging q 
and p. 

We will consider classical symbols a £ S®(1). As was stated in Section 
2, Weyl quantization of functions in Sq(1) leads to bounded operators. The 
following result sharpens this statement. 

Proposition 5. [Calderon-Vaillancourt] There is a C < oo and a finite n £ IN 
such that for all a £ Sg(l) and £ € [0, 1] 



|a|<ra,|,3|<n,(a;,p)GlR d xM 



(30) 



For a proof see 0, Theorem 7.11. The statement there is slightly different, 
but their proof implies our Proposition]^. Note, in particular, that the constants 
n and C in Proposition [5] depend on the dimension d of configuration space. 

The so called product rule, presented in the following Proposition]^, is at the 
basis of all semi-classical analysis we will apply. Given Proposition |3|, however, 
its proof is mainly a matter of calculation. For details see |||, Proposition 7.7 
and Theorem 7.9. 

Proposition 6. [Product rule] Let a,b € So(l). Then for all n £ INo there is 
a d n < oo such that 



a W,e 6 W,e_ ( £ £fccfc) 



W,e 



£(i 2 (M)) 



< d n e 



n+l 



(31) 



with 



Ck{x,p) 



E 

|a| + |/3|=fe 



(_1^ 



({dga*a)(agd5b)){x,p). (32) 



We state two simple facts that are immediate consequences of the product 
rule as 



Corollary 7. Let a, b G 5g(l). Then 
(i) 



0( £ ). 



£(L2(M)) 

(ii) If, in addition, supp(a) n supp(fe) = 0, then for any n £ IN 



a W,e b W,e 



£(L 2 (M)) 



0( £ ") . 
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Proof. (i) is just Proposition |6| for n = 0. (ii) holds since in this case = 
in @ for all fc <E 1N . 

□ 

Here and in the following 0(e n ) means that an expression, or its appropriate 
norm, is bounded by a constant times e n for sufficiently small e. 

The crucial ingredient to our semi-classical analysis is the following first- 



order version of a Theorem going back to Egorov |3j (cf. also 12 , Theoreme 
IV-10), which is also a direct consequence of the product rule. 

Proposition 8. [Egorov's Theorem] Let a g Sq(1) and < T < oo. There is 
a C < oo such that for all t S [-T, T] 



Ui 



Ul 



(00$') 



<Ce 2 



£(L 2 (M)) 



Proof. First note that for a € S^ 1 ) we have that a o $' 6 Sq(1) for all 
finite i, since the Hamiltonian vector field is smooth, uniformly bounded and 
all its partial derivatives are uniformly bounded. Moreover, a o $ and all 
its partial derivatives are each bounded uniformly for t € [—T,T]. Therefore 
|j (a o $*) W ' e ||£(L2( M )) is bounded uniformly for t £ [—T,T] by Proposition |[ 
Writing 



t-s\W,e 



one is led to consider 



(33) 



(34) 



r,(- 

W,E 



where {•, •} denotes the Poisson bracket and we used that Hi = H , ,£ . From 
the product rule one computes that for arbitrary a, b G Sq(1) 



)={a,6} W < £ + 0( £ 2 ), 



(35) 



which implies that (|4|) is 0(e 2 ) for fixed t — s. However, one can easily convince 
oneself that this bound is uniform for t — s in any bounded interval, since 
the semi-norm used in Proposition |5| is bounded uniformly for the C2-terms 
appearing in a derivation of Equation ( |35| ) using the product rule. 

□ 

To have a natural way for extending functions of x or p alone to functions 
on phase space, we introduce the projections tt x : H d x M — > JR d and tt p : 
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H d x M — > M as tt x {x,p) — x and ■K p (x,p) = p. We are now ready to prove 
Lemma [lj. 

Proof (of Lemma II]). In order to regularize t Ai and we pick some f\ G 



C^{M) and f m G C°°(M) such that fA 



1, /n 



I Af\A„ 



1 and supp(/i 



$ p *)n supp(/ m o 7T P ) = for all t £ [0,T^]. Egorov's Theorem implies that 
there is a C < oo such that for t G [0, ZJ°J 



(/i°7r p ) w ' e f/i 



(/i ° 



t\W,e 



<Ce 2 



C(L 2 



Since, by construction, supp(/i o $ p *)n supp(/ m o ir p ) 

-t\W,e 



5, we have (f m o 



7r p ) w ' e (/i o $-*) W ' e = 0(e 2 ) by Corollary g. This and the fact that P; 
(/i o 7r p ) w - £ Pi and that (1 - P m ) = (1 - P m )(/ m o vr p ) w ^ implies (§" 



□ 



4 Convergence of the unitary groups 

We prove Theorem |[ Let A; + S C A g . In the following T < oo with T < 
T^(Ai,A g ) will be fixed once and for all and we will always assume that < 
t < T. 

For reasons that will become clear during the proof we have to introduce 
further sets in momentum space. Let 

A m = (J supp(]l Ai o$ p *(>,•)) 

t£[0,T],xeH d 

be the subset in momentum space that is reached by the classical dynamics 
starting from A; in our relevant time span. By construction we have that A m + 

5 C A g and we define A mi := A m + <5/4, A m2 := A m +8/2 and A m3 := A m + 35/4. 
The respective projections will be denoted by P m i, P m 2 and P m 3. The following 
Lemma ensures that V s maps wave functions supported in A m i (resp. in A m 2, 
A m 3) to wave functions supported in A m 2 (resp. in A m 3, A g ) up to errors of 
arbitrary order in e. We will use this result in the following implicitly many 
times. 

Lemma 9. Let V G L 1 ^) such that / dk \k\ n \V(k)\ < oo for all n G IN. 
Then for each compact A C K d and all m G IN and 8 > there is a C$ m such 
that 

UV S -P A+S V^P A \\ C{L2{M)) < C s , m s m . 

Proof. Let (f> £ RanPA, i-e. supp</> C A, and note that 
{V £ 4>){p) = I ' dkV{k)<j){p-ek) 
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dkV{k)(j){ P -ek)+ J dkV{k)<f){ P - ek) . 

fc|<e-!/2 \k\>e~^/ 2 

(36) 



The first term in ( |36[) is supported in A + S for e sufficiently small. For the 
second term note that 

dk V(k)(f>(p -ek)= J dk' e- d V(k'/e)<j)(p - k') 

amounts to convolution with the function l| fe ,i >e i/2£ je) and therefore, 

by Young's inequality, the £(L 2 )-norm of the corresponding map is bounded by 

dk'e- d \V(k'/e)\= J dk\k\- 2m \k\ 2m \V(k)\ < C m e m . 

W^e 1 ' 2 \k\ye~ 1 ' 2 

□ 

We turn to the proof of (|23|). Using Lemma [j] together with UU* =Id on 
RanP m 2 we get 

* (;)-«•* (;)")«- 

= ~ iU f~) Jj dsU{-s) {HU* - U*Hi)Ui{s)U Pi 

= -iU J° dsU(-s) (H - P m2 U*Ui(s)U Pi + 0(e) . 

(37) 

For the last equality note also that a factor of order 0(e 2 ) in the otherwise 
uniformly bounded integrand leads to the integral being 0(e). We would be 
done by the same argument, if we could show that H —U* H\IA acting on RanP m 2 
is 0(e 2 ). However, the first order term does not vanish, as we will see, and we 
have to treat the integral more carefully. 

In order to separate the leading order term we write 

H - U*H X U = {H- H dias ) + (ff diag - WH{U) , 

where 

H diag = H + P g V £ P g + P^V e P^ (38) 

and P^ := 1 H - P g . 

We will treat the easy part first and show in Lemma [n] that the difference 
-ffdiag —U*HiU vanishes sufficiently fast on RanP m2 - 
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However, to do that we need to show that Po(p) is twice continuously diffcr- 
entiable as a function of p and, for later purposes, we give an explicit expression 
for V p P (;p)- 

Lemma 10. P (-) £ C 2 (A g , C(Hf)) and for p £ A g one has 

(V P P )(P) = -(ReQo)(p)(V p H )(p)P ( P ) - P (p)(X7 p H )(p)(R E Q )(p) , 

(39) 

where {R E Qo)(p) = (#o(p) - E{p)) X Q Q {p) and Q„(p) = l Ht - P (p). 

Note that (ReQo)(p) is bounded since £?(p) is an isolated eigenvalue for all 
p £ A g and Qo{p) projects on the orthogonal complement of the corresponding 
eigenvector. 

Proof. For p £ A g we can express Po(p) as a contour integral: 



Po(p) = ~{ dXR x (H (p)), 

lm Jc(p) 

where c(p) is a smooth curve in the complex plane circling the isolated eigenvalue 
E{p) only and R x (H (p)) = {H (p) - A)" 1 . Multiplying 

= V p l = V P (H Q (p) - X)R x (Ho(p)) 

= (V p H (p))Rx(H (p)) + (H Q (p) - X)(X7 p R x (Ho(p))) 

with R x from the left we get 

\7 p R x (H (p)) = -R x (H (p))(V p H (p))R x (H (p)), 

which is bounded uniformly for A £ c{p) by our assumptions on the family 
Hq(p). Hence 

V P P ( P ) = I d\V p R x (H (p)) . 

Jc(p) 

Collecting the above information we obtain 

Qo(p)V p P (p) = Qo(p)V P Po(p)(Po(p) + Qo(p)) 

= -LI d\Qo(p)R x (H (p))(V p H (p))R x (H (p))P (p) 

lm Jc(p) 

= -Re(H (p))Qo(p) (V p ff (p)) Po(p) , 

where the Qo(p)^ p Po(p)Qo(p) vanishes, since in this case the integrand is an 
analytic function inside of c(p). Therefore Qo{p)^ p Po{p) = Qo{p)^ P Po(p)Po(p) 
and, using self adjointness of V p P (j>), it follows that 

Po(p)V p P (p) = Po(p)V p Po(p)Qo(p) = (-Re(H (p))Qo(p) (V p H (p)) P (p))* . 
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This yields ( |39| ) and Pq(-) G C 2 (A g , C(Ht)) by applying the same argument 
again to (|39|). 

□ 



Lemma 11. The phase of tpo(p), P 6 A g , can be chosen such that V'o(') € 
C 2 (A g ,7if). Then for e sufficiently small there is a C < oo such that 



[H^-WH^lAPva <Ce 



(40) 



Proof. Note that Lemma |^ implies that 

\\PgH di , g P ni2 \\ < Ce\ 

and we have by construction that P g U* HiUP m 2 — 0. Hence it suffices to 
consider the difference projected onto RanP g . On RanP m 2 we have 

P s (H dl3 . s H> )(p) = 

E(p)(b( P )Mp) + U e (p) J dk9(k)<i>(p-Ek) (Mp),Mp-ek)) Ht Mp) 

and 

{WH x UH> ){p) = E(j>)<t>(p)i> (j>) + 1a, (p) / dk V(k)cf>(p - ek) Mp) ■ 
Hence 

P g (i?dia g - U*H X U) (<M>)Cp) = (41) 
= 1a, (p) f dkV(k)<P(p-ek) «Vb(p),Vo(p-e*))w t -l) V>o(p) ■ 
We will show that there is a constant C such that 

K^)(p),Vto(p-eA))w t -l| <C\k\ 2 e 2 (42) 
for all p € A g and k with p — ek £ A g . Therefore 

(H diag -U* 11^)^0^ < Ce 2 J dk\V(k)\\k\ 2 \\cf > (--ek)M-)\\ H 

< Ce 2 U\\ L 2 {M) J dk\V(k)\\k\ 2 

= Ce 2 \\Mn [ dk\V(k)\\k\ 2 . 



To see note that for ^o(') G C 2 (A g ,Hf) Taylor expansion yields 



^o(p - efc) = i>o(p) - ek ■ V p Vob) + \e 2 k ■ H(ip )(jp')k , 



17 



where H{\p ) denotes the Hessian and \e 2 k ■ H(ip )(p')k is the Lagrangian re- 
mainder. In view of (tpo(p), V F ^o(p))H f — 0, which follows from comparing ( |39| ) 
with 

(V p P ip){p) = (tpoip), V>)WfVpVo(p) + (Vpipoip), ip)n t ipo(p) , 

we obtain 

|<Vo(p), ^(p - £fc))« f - 1| < ^b)|fc| 2 £ 2 ■ 

Here C(p) = | £V j KV'o(p'), c^c^ V'oCp'))!; which is bounded uniformly for p £ 
A g by continuity. 

We are left to show that one can indeed choose the phase of the eigen- 
functions such that tp (-) € C 2 (A g ,7Y f ). Since P (-) <E C 2 (A g , C(H { )), one 
can cover A g with finitely many open sets U n , each containing a p n such that 
ipoip) = Po(p)ipo(Pn) /\\Po(p)i>o(P7i)\\ is well defined for all p E U n . It is now 
straightforward to connect these pieces in order to get a C 2 version of V'o(p) on 
all of A g . Note that if we would replace smoothness by analyticity, this last step 
would become nontrivial, in particular, if M is a torus (cf., e.g., ]l0(]). 

□ 

With the help of Lemma [ll] and ([37j), we have 

= -iU J ' ds U(-s) {H - P diag ) P m2 WU X [s)U Pi + 0(e) . 

In the following lemma we isolate the leading order term in (H — i?<jiag)Pm2- 
Lemma 12. 

(H - P diag )P m2 = -eQ g (V p P g )P m3 ■ P £ P m2 + 0(e 2 ) , (43) 
where (V p P g ) = dpV p P (p), Q g = dpQ (p) and 

(P £ V)(p) :=-« / dkkV{k)ijj(p-ek). 



Proof. Using Lemma g we obtain 

(P - Pdiag)Pm2^ = Pg^P^ 

rfp ^Qo(p) y dfc F(fc)P (p - efc)(P m2 V>)tp ~ e*0 ) - < 



Since Pq : A g — > C(Tii) is twice continuously differentiable, we have that 

P (p - efc) = Pob) - ek • (V p Po)(;p) + e 2 k ■ H(P )(p'(p, ek)) ■ k , (44) 
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where the last term is the Lagrangian remainder with H denoting the Hessian. 
Hence, for p e A g , 



dk V(k)P (p - ek)(P m2l p){p - ek) 

dkV(k)(P (p) - ek • (V p P o )(p))(fm2V0(p ~ ek) (45) 
+ e 2 [ dk V(k) k ■ H(P )(p'(p, ek)) ■ k (P m2 V)(p - ek) . (46) 



Since 



lA m3 (-) / dkV(k)k- H(P )(p'(;sk)) ■ fc(P ro2 V)(- - ek) 



(47) 



H 



< J dkV(k)\\i Ani3 {-)k-H(P ){p'(-,ek))-k{P m2 ^)(-~ek)\\ n (48) 

< sup \\H(P )(p)\\ f dk \V(k)\ \k\ 2 \\l Am3 (-) (P m2 V)0 - Ek)\\ n 
peA m2 J 

< cu\\ H I dk\v(k)\k 2 , 



(46) is 0(e 2 ) in C(H{) and multiplying ( [45| ) with Qo(p) from the left establishes 

□ 

Including Lemma [l^ we are left with 

= leU £ dsU(-s)Q g (V p P g )P m3 • F e P m2 Zrf/i(s)WPi + 0(e) . 

To exploit the time averaging we write <5 g (V p P g ) as a time derivative, at 
least in approximation. Let 

B( P ) := RE (p ){H a {p)) 2 Qv(p){V p H a )(p)P a {p) . 



Lemma 13. 



and 

in £(H). 



5 (V p P g )P m3 = [B,H ]P m3 

= [B,H]P m3 + 0(e) 
[P £ ,P]P m2 =0(e) 



(49) 
(50) 

(51) 
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Proof. Recalling (|39j), we have for p £ A g that 

Qo(p)(V p P )(p)Po( P ) = -R E(p) (H (p))Q (p)(V p H Q )(p)P (p) 

and thus ( fl^ ) follows from direct computation, 

B(p)H (p) - H a (p)B(p) = 

= -(Ho(p) - E(p))R E{p) (H Q (p)) 2 Q (p)(V p H )(p)P Q (p) 
= -R E{p) (H (p))Q Q (p) (V p H Q )(p)P Q (p) 
= Qo(p)(V P Po)(p)Po(p). 

For @ we need to show that [B, V £ ]P m3 = 0(e): 

[B,V s ]P m3 = [R E (H ) 2 Q s (V p H )P g ,V e ]P m3 . 

It follows from the proof of Lemma |l2| and using Lemma |] that 

[P g , V e ]P m3 = P g V £ P m3 - P g V £ P g P m3 - Q g V £ P g P m3 + 0(e) 
= -Q g V £ P g P m3 + 0(e) = 0(e) . 

Hence 

[B, V £ ]P m3 = [R E (H ) 2 Q g (V p H ], V £ ]P g P m3 + 0(e) . 

Using again Lemma ^| and the fact that also W p H acts on p-fibers, we observe 
that 

([V p H Q ,V £ ]P m3 ^)(p) = t Ag (p)([V p H ,V s }P m3 iP)(p) + O(e) 
= t Ag (p) J dkV(k)((V p H )(p)-(V P H )( P -ek))(P m3 i;)( P -ek) + O(e) 

= et Ag (p) [ dkV(k)k-H(H )(p'(k))(P m3 ^)(p-ek) + O(e) (52) 



where H(Hq)(p' (k)) denotes again the Hessian evaluated at the appropriate 
point p'(k). By assumption we have that H(Hq)(p) is bounded uniformly for 
p £ A g in the graph norm, which is equivalent to the 7i-norm on span{?/;o(p) , p € 
A g }. Thus the above expression is 0(e) as e — * (cf. estimate (0)). Now 

[B, V £ ]P m3 = [R E (H Q ) 2 Q g , V £ ](V p H Q )P m3 + 0(e) . 

Since R E ( p )(Pto(p)) 2 Qo(p) is bounded and differentiable with respect to p, one 
can show by the same line of arguments as in the case of [P g , V £ ]P m 2 (cf. proof 
of Lemma |l2|) that 

[R E (H ) 2 Q g ,V £ ](V p H )P m3 = 0(e) . 



Thus we showed @. For (|I|) first observe that [F e , V £ ] = 0. [F e , H ] = 0(e) 
follows from an estimate analogous to ([52]). 

□ 
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Let A — B ■ P m3 F £ . Then Lemma [jj yields 

Qg(V p P g ) • P m3 F £ P m2 = [A, H]P m2 + 0(e) 

in £(ft), where [P m3 F e } H}P 2 = 0(e) follows immediately from @ Let A(t) 
U(-t)AU(t), then 

= ieU (^j dsU(-s)[A,H]U(s)U(-s)P m2 U*Ui(s)U Pi + 0(e) 
= -eu(-) f S ds (-^-A(s))u(-s)P m2 U*Ux(s)UP i + 0(e) 



-eU\- 



_-r | - 



e / ./o 

A(s)U(-s)P m2 U*Ui(s)U ° Pi + 0(e) 



(i 



ieU 



dsA(s) — (u(-s)P m2 U*U 1 (s)u) Pi 
° dsA(s)U{~s)(HP m2 U* - P m2 U*H^j U x (s)UPi + 0(e) 

jf 6 ds A(s)t/(-s) (i? - WHiUj P m i U*Ui(s)U P ; + O(e) . 

For the last equality we used again Lemma |l| and the fact that Lemma ^| guar- 
antees P m2 U* HiUP m \ = U*HiUP m i + 0(e 2 ). Finally also the last integral is 
0(e), since we showed already that (H — U* HiU)P m \ — 0(e) and the proof of 
Theorem || is completed. 

5 Convergence of the macroscopic observables 

In this section we prove Theorem [| and Theorem ||. We start by showing 

Lemma 14. Let g e C°°(M) such that goir p E Sq(1) and let the assumptions 
of Theorem ^ be satisfied. Then 

\(u(-t/e)(g(p)®l)U(t/e)-U*U 1 (-t/e)g(p)U l (t/e)U)p\< Ce, (53) 

where g(p) := (go7T p ) W ' £ denotes the operator of multiplication with the function 
g on L 2 (M). 

Proof. Let A i>7 = Aj + 7 and P i>7 the corresponding projection. Let 7 be 
sufficiently small to ensure that supp((lA ij7 o •)) + 6/2 C A g for all 

t G [0,T*J and x G R d . We abbreviate g(i) := U(-t/e)(g(p) ® l)U(t/e) and 
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gi(t) := Ui(—t/e)g(p)Ui(t/e) and split J53|) into two parts 



(g^-Wg^U^P, 



< 



P^Jg^-Wg^u)^ 



(54) 
(55) 



Lcmma| 



< 



Using UU*P g = P g on L 2 (M) we find that Q) can be estimated as 
P in (g{t)-W gi {t)u)P^ 

Pi n (u {-t/e)U* - U*V X (-t/e) ^UU*P g g(p) U x (t/e) UPi 

+ jp ii7 U(-t/e) (ff(p)®l) {u{t/e)-WU 1 {t/e)U)P 1 
Pi,y(u (-t/e) - WUi (-t/e) w) P g g(p) U t (t/e) UP, 



0(e 2 



0(e 2 



Thcorc 



+ 

0(e 



P i)7 t/(-t/e) (5(p)®l)P g | \(u(t/e)-WU 1 (t/e)u)P i 



It remains to show that in ([35]) both terms are of order 0(e) separately. By 



I Aj 



1 

' Jin 



7 I M\A ii - I 



construction we can pick /i, /j-, e C°°(M) such that /i| 
and supp/iP| supp/^ = 0. Using Egorov's Theorem and Corollary [?] (i) we get 

5iW(/i°T P ) W ' £ = ( 5 °^) W ' e (/iO^) W ' E + 0( £ 2 ) 
= (07o$*)(/ 1 o % ))W.« + O(e ) 

and therefore, with Corollary [7] (ii), that 



\W,S ; 



||(/iW) w ' £ (G?°$p)(/i°%)) 
0(e), 



U Pill + 0(e) 



since f^- o ir p and (g o $* ) (/; o n p ) are disjointly supported. Finally we compute 
that 



Pdr9(t)Pi\\ 


= \\ P ^ U (- 


-t/e) {g(p)®l)U (t/e) 






Theorem [j] 


\\ P ^ U (- 


-t/e) U* g(p)Ui (t/e)UPi\\ + 0( 


s) 




Egorov 


\\ P tf U (- 


-t/e)WU 1 (t/e) (jo^pf/io 


7T p ) W ^UPi 


|+0(e 


Corollary ^J(iJ 


\\ P in U (- 


-t/e)WU 1 (t/e) (/i°%) W ' E (<?° 


$l) w ' e UPi 


|+0(e 


< 


\\ P uyU(- 


-t/e)U*U l (t/e)UP in U*(go$l 


) w ' e UP { \\ + 


0(e) 


Theorem [j] 


\\ P ^ U (- 


-t/e)U(t/e)P- in (go^l) w ^UP i 


I + 0(e) 






0(e). 
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For the third equality we inserted Ui(t/e)Ux(—t/e) = 1 in order to apply 
Egorov's Theorem on g\ (t) . 

□ 



We now come to the proof of . We have (at the moment only formally) 



that 



x E (t) = ieV p <g> 1 + / dsU(-s/e)[iV p ®l,H]U(s/e) 



On HiQR?) ® W f this gives 

x E (t)Pi = (ieV p ®l)Pi+ / dsU(-s/e)[iV p ®l,H}U(s/£)Pi 

Jo 



(ieV p <e>l)Pi + / 
Jo 



ds U(-s/e) (V p #o) P g ^(s/£) P 



+ / dsU(-s/e)(V p H )P^U(s/e)P i 
Jo 



= (ieV p ®l)P+ / dsU(-s/s) (iVE(p)®l)P s U(s/e)P 1 
Jo 

+ 0(e), (56) 

where we prove in a moment that this is well defined on H 1 (JR d ) (S>?if and that 
the last equality in ( p6] ) holds. 

On the other hand, again on H 1 (M d ) £g> 7if, 

U* x\(t)UP i =U* ieVpUPi + [ dsU* Ui(-s/s)iVE(p)U 1 {s/e)U P. 

Jo 

Since VP € C yD (M) we can apply Lemma [l4| and finish by showing that 

(ieVp® 1-WieVpUjPi = 0(e) . (57) 
Let ip = (f>(p)ipo{p), then 



(i eV p ® 1 - It i eV p Z-/) 0Vo 



H 



< eC||0|| L 2 (M) =eC\\ip\\-H ■ 



which proves ( |57|). To complete the proof of Theorem || we are left to show the 
last equality of (|56D, i.e. that 

(V p P )P g i C/( S /£)P i = O( £ ) 

uniformly for s <E [0, T\. If (V p Po) was a bounded operator, Lemma |l| translated 
to the full dynamics via Theorem ^ would take care of that. Recall, however, 
that by assumption (V p Po)(Po — *) _1 is bounded and hence it suffices to show 

H Q P^U(s/e)P i = 0{e). 
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We calculate 



H Q P^U{s/e) Pi = P£ H U{s/e) P x 
(*) 



Lemma 9 
Lemma 







P£ HU(s/e)P i + 0(e) 
P^U(s/e)HP i + 0(e) 
P^U(s/e)P in HP i + 0(e) 
0(e) , (58) 



where (*) follows from Lemma [j] together with [P^~, V e ]P m = 0(e) (cf. proof of 
Lemma [D]) and P ltl is defined as in the proof of Lemma [l4| 

Finally we sketch the proof of Theorem [I| First we conclude from Theorem 
U that in the special case of a = / o tt x G S®(1) H Coo(]R d x M), for some 
/ € C°°(TR d ), ( |28[ ) holds. To see this, note that by the functional calculus for 
self-adjoint operators it suffices to show 

\\(f(x s (t))-f(Wxt(t)U))P i \\<Ce. (59) 

However, ( |59| ) follows from a standard approximation argument using the gen- 
eral Stone- Weierstrafi Theorem and norm-convergence of the resolvents (cf., e.g., 
Theorem VIII.20 in fujl). 

Next we observe that for a(x,p) = f(x)g(p) Lemma [j] and Corollary ^ imply 

in a rather straightforward way. 
Finally one can approximate general a G Sq(1) H C 00 (lR <i x M) by products 
referring again to the general Stone- Weierstrafi Theorem. 



6 Semi-classical distributions 

Theorem H establishes that, restricted to the ground state band, the full unitary 
group U(t) and the approximate one-particle unitary Ui(t) are uniformly close to 
each other and Theorems || and ^ lift this assertion to semi-classical observables. 
Experimentally measured are empirical statistics of suitable observables, like 
position and momentum, and we still have the task to investigate in what sense 
they are approximated through the time evolution U±(t) for e< 1. 

For small e, H\ itself is a semi-classical Hamiltonian, which means that em- 
pirical distributions can be determined through the classical flow $ generated 
by H\. Somewhat crudely the scheme is as follows: one choses an initial wave 
function ip E , which may or may not depend on e, such that for small e it deter- 
mines the measure p c \(dxdp) on phase space. We evolve ip £ as = e~ tHlt l e iJf 
and p c \(dx dp,t) = (p c \ o $>~ t )(dx dp). Then the empirical distributions com- 
puted from ipt agree with those of p c \(t) up to errors of order e, i.e. quantum 
distributions are well approximated by their classical counterpart. In our con- 
text the true statistics must be compared with U(t/e)ip e and we have to make 
sure that the approximations of Theorem [| are so sharp that we can draw the 
same conclusions for U(t/e)ip £ . 
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There are various "schools" which differ in what initial ip's are regarded as 
physically natural. 

(i) wave packet dynamics. The initial wave function is well localized in macro- 
scopic position and momentum, i.e. p c \(dxdp) = 8(x — Xo)S(p — po) dx dp. Then 
the wave packet follows the classical orbit , which only reflects that p c \(dxdp, t) = 
S(x — x t )S(p — pt) dx dp. 

(ii) microscopic wave function independent of e. On the macrscopic scale the 
position is localized, but there is momentum spread. Therefore p c \(dxdp) = 
8{x)\^{p)\ 2 dxdp. Such a choice is appropriate immediately after a scattering 
event. Then ip is still localized at the scatterer but has considerable momentum 
spread. 

(iii) WKB. The wave function is taken to be build up from local plane waves, 
which means it has the form ij) £ (x) = e d l 2 f(£x)e lSty£X ^ e on the microscopic 
scale, ip 6 is spread over macroscopic distance, but at any given point it has a 
sharp momentum. ip £ yields the phase space measure p c \{dxdp) — \f(x)\ 2 S(p — 
VS{x))dxdp. 

As will be discussed, all three classes of initial wave functions, and in princi- 
ple more, can be handled under rather mild regularity assumptions in a unified 
fashion. The approximation from U(t/e) to U\{t/e) is covered by Theorems || 
and [| and the semi-classical distributions for U\ (t) are a consequence of stan- 
dard results of semi-classical analysis as presented, for example, in [Q. Since we 
did not find a discussion of sufficient generality in the literature, we explain the 
arguments in some detail. 

Note that in the following we only treat the case M = M d . If M is a torus, 
position and (quasi)-momentum are not related by Fourier transformation, but 
by a Bloch-Floquet transformation. Apart from this difference the analysis goes 
through analogously. 

6.1 Wave packets following a classical trajectory 

The conceptually simplest way for a quantum particle to behave classically is 
to have a well localized wave function that follows a classical trajectory. Hence 
we consider initial wave functions with sharply peaked momentum and, at the 
same time, sharply peaked macroscopic position. Let the initial wave function 
be 

for some <f> € L 2 (fR d ), i.e., a wave function that is peaked on the macroscopic 
scale and centered at Xq, but spread out on the microscopic scale. Its Fourier 
transform is given by 

2 I \ -d/4 -i£oi£^£oj 2tP ~P0\ 

0x o ,po(p) = £ ' e E <P( — T-), 

V e 

which becomes sharply peaked at po for e small. There is no difficulty to include 
also asymmetric scaling with weights e 1_Q and e a , < a < 1, and the choice 
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a = 1/2 was made to simplify presentation. Under the time evolution generated 

by 

Hi =E{-iV x ) + V{ex) 

it moves along the corresponding classical trajectory starting at (xojPo) follow- 
ing the classical flow $* generated by 

H cl = E(p) + V{x) . 

To be consistent with the previous chapters, we continue to work in momentum 
representation. 

Proposition 15. Let a E % G ff 1 (H d )nD(p) and T < oo. Then there 

is a C < oo such that for t G [0, T] 

< cVi||?||(||0|| + WI + ||v^||). 

(60) 

Using Theorem ^, this translates immediately to the full dynamics. 

Corollary 16. Let the assumptions of Theorem |] be satisfied. Then there is 
a C < oo such that for ip £ ((H 1 ^) n D(p)) ® W f )) n RanP; and t € [0, T] 

| (ipx , P0 ,a E (t) ip Xo>Po ) - (ao^ t )(x ,po) | 

< cyi||v>||(M + ||(p® i)^|| + 1| (v p ® i)vil) ■ 

Hence, when initially, on the macroscopic scale, position and momentum 
are both sharply defined, the wave packet follows the classical orbit without 
spreading even for macroscopic times. Such a situation occurs for example in 
particle accelerators, where one can indeed calculate the particle trajectories 
based solely on classical dynamics in good approximation. 

Proof [of Proposition ^5[] Referring to Egorov's theorem we have to compute 

^, p „(ao$ f ) w ^;, pt ). 

One can replace (a o $*) w > £ by the so called standard quantization of a o $* 
defined by 

((ao$*) s - e 0)(p) := (27r)- d/2 J di(flo$ f )(a, P )e- ,p ^(i) (61) 

where the error is of order e uniformly in t G [0,T] (cf. Chapter 7 in B). For 
the standard quantization the result becomes a simple calculation: 

(?zo,Po>(a° $i ) S ' e ?*o,Po) = 



(?zo,po> a l( t )</Wo) - ( ao **)0°'Po) 
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= (27r)- d / 2 I dxdpe*"^* 1 ^ (jjr) ( ao ® t )( ex >P) e ~ W ' Xe * X ' P ° x 

= (2tt)-^ 2 J dxdp^* (^j {ao^ t )(ex + x ,p + P o)e~ l P- x ^(V^{x)) 
= (2n)- d/2 J dxdp^*(p)(ao^ t )(y^x + x ,^p + po)e- t ^(t)(x) 

= (ao$*)(a; ,p ) + O(VI), 

where the last equality follows from a Taylor expansion of ao$ f around (xo,po) 
together with the assumption that <f> <E H 1 (fii d ) n D(p). Note also that (V K (ao 
$*)) £ Sg(l) and (V p (ao $')) e S°(l) and hence (V x (a o $*)) S > E , (V p (o o 
<E>*)) S < £ € £(L 2 (IR d )) (cf. Chapter 7 in §) 

□ 



6.2 Initial wave function with momentum spread 

Generally wave functions are not of the special form described in the previous 
subsection. Typically they arise from microscopic interactions and thus "live" 
on the microscopic scale and do not depend on e. But if the shape of the 
initial wave function does not depend on e it will effectively look like a delta 
function on the macroscopic scale at t = 0. More precisely, let 4> € L 2 (fR d ), 
then the scaled position distribution is e~ d \4>(x/ e)\ 2 which converges to S(x) as 
a measure. However, the scaled momentum distribution is still |</>(p)| 2 , since in 
the quotient x/t the e's cancel and thus the initially peaked wave function will 
spread if evolved to times of order e~ l . More generally we consider as initial 
wave function 

<Px {x) = <j>{x - x /s) , 

i.e. we move <p to the macroscopic initial position xq. 
Then it is natural to chose 

p c i(dxdp) = 5{x - x )\(f) Xa (p)\ 2 dxdp 

as the corresponding classical phase space distribution at t = and evolve it 
according to the classical flow to 

p c \{dx dp, t) = (p o <j>~') (dx dp) . 

Let us now, as the simplest example, compare the quantum mechanical po- 
sition distribution 

Pe {dx,t)=e- d \<t> t {x/e)\ 2 dx, <j> t = e- lH ^l^ XQ , 
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with the classical one, 

pl\{dx, t) = J pd(dx dp, t) , 

in the limit e — > 0. As a first step we calculate, with / £ Cfi° a test function, 

J p e (dx,t)f(x) = {fi,f<%) = (4> xo ,f(x s (t))<p xo ) 

= <0x o ,(/o$*) W ' £ ^o)+O( £ 2 ). (62) 

For the last equality we used Egorov's Theorem. 

To proceed we need to know how the Weyl quantization of a time evolved 
classical observable acts on microscopic wave functions. We postpone the proof 
of the following Proposition to the end of the section. 

Proposition 17. Let a £ Sq(1). Then for each T < oo there is a C < oo 
such that for 4> £ D(x) and t £ [0, T] 

< Ce\\xct)\\ , (63) 

where (a o $*j(xo, •) denotes the operator of multiplication with the function 
(a o $*J(xo,p) in momentum representation. 

Thus ( |62| ) becomes 

(?xo,(/°*x) W,e ?xo>= I dp\d X0 (p)\ 2 (fo^ x )(x ,p) + O(s). (64) 
However, the right hand side of (|64| ) is exactly what we were aiming for: 
J dp\^ xo (j>)\ 2 {fo^ x )(xo,p) = J p d (dxdp,0)(fo^ x )(x,p) 

= J p c \(dxdp,t)f(x) 
= J pUdx,t)f(x). 

In summary, one obtains that 

lim p e (dx,i) = Pci(dx,t) (65) 

e— *0 

weakly as measures. Note that this means that the position distribution of the 
quantum particle converges to the classical one, although the wave function 
is not following the classical orbit, but spreading. Such a situation occurs for 
example in a scattering experiment. After the quantum particle is scattered off 
the target its wave function usually has a large momentum spread and thus it 



((«o$f £ -( tt o$')( Io .))^ 
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also spreads in position space. However, during the process of detection it is 
subject to relatively weak potentials and can be treated like a classical particle, 
at least on the level of statistics. 

Such a result holds in fact not only for the position distribution but for all 
semi-classical observables. For general observables a G Sq (1) we get analogously 

{<f>x ,al(t)<f> Xo } = J p c i(dxdp,t)a(x,p) + 0(e) . (66) 

Now Theorem [I] allows us to immediately translate this result to the full 
quantum dynamics. 

Corollary 18. Let the assumptions of Theorem || be satisfied. Then there is 
a C < oo such that for ip G [H 1 (H d ) <g> Hi) n RanP ; and t G [0, T] 

< CeW(W + ||(V p ®l)V»||). 



(ip Xo) a s (t)ilj Xo ) - / p cl (dxdp,t)a(x,p) 



Before we close the section with the proof of Lemma LL7L let us shortly com- 
ment on an approach "dual" to Weyl quantization of classical observables. One 
can translate the above result to the level of distributions, i.e. generalize ([35]), by 
looking at the scaled Wigner distribution. For <fr £ L 2 (fR d ) it is defined through 

W £ (dxdp,t) = J die- d (t)* Xo {x/e-i/2,t) <t> Xo {x/e + f/2,t) e ip< dxdp. 

The scaled Wigner distribution yields expectations of Weyl quantized operators 
through 

(tj> xa , a\(t)4> Xo ) = J W £ (dxdp,t)a(x,p) . 
This together with (|6^) implies 

lira W £ (dx dp, t) = p c \{dx dp, t) 

s— >0 

weakly as measures. 

On the level of the full quantum dynamics we introduce the reduced Wigner 
distribution for ip G L 2 (Si d ) ® Tit as 

W £ d (dxdp,t) = J dte- d {r xo {x/e-i/2,t),^ Xo {x/e + d/2,t)) H[ e^dxdp, 

where here ^ stands for Fourier transformation in the first argument, i.e. for 
T ® 1. Corollary [l8| translates into 

lim W £ d (dx dp, t) = p c \(dx dp, t) 

weakly as measures. 
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Proof [of Proposition [l7|] . Again we replace (ao$ f ) W ' c by the standard quan- 
tization ofao$ f and get 



((«o§') s ^)(p) := {2n)- d ' 2 J dx(ao^) (ex,p) e*>* 4>(x - x /e) 

= (2n)- d ^ [ dx (a o $*) {ex + x ,p) e~ l ^ x e -*f *°/ e ^ . ( 67 ) 



Taylor expansion yields 

(oo $*)(ea; + .To,p) = (ao$')(i„,p)+£i • (V s (a o $ t ))(y(a;),p) . 



We insert this into (|6j]) and conclude as in the proof of Proposition |15| that the 
term proportional to e is bounded in norm by our assumption that <f> £ D(x). 
This proves (p3|). 

□ 



6.3 Initial wave function of WKB form 

In the previous case of an initially localized wave function the different momen- 
tum components travel at different velocities and therefore such a wave function 
spreads on the macroscopic scale. After some time one expects the wave func- 
tion to have locally well defined momentum as long as no interference occurs, 
i.e. it should be of WKB type. On the microscopic scale a WKB wave function 
has the form 

<j>{x) = e d ' 2 f{exy— , 

with / and S real valued. Hence it locally looks like a plane wave with momen- 
tum VS(ex) and amplitude f(ex). 

Time-dependent WKB approximation is concerned with showing that the 
time evolution of such a wave function is in first order given by 

where St is the solution of the classical Hamilton- Jacobi equation with initial 
condition S and ft is the solution of the classical continuity equation 

d t f + divVSt = 0. 

The corresponding classical phase space distribution is therefore 

p c i(dxdp) = f 2 (x)S{p-\7S(x))dxdp. (68) 

The main drawback of the WKB approximation is that it works only as long as 
no caustics are reached, or, put differently, as long as no interference between 
different parts of the WKB wave function happens. 

However, if we focus again on distributions of semi-classical observables on 
phase space no difficulty arises, only p c i(t) is no longer of the particular form 
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Proposition 19. Let / E C°°{1R d ) n L 1 (]R d ), S £ C°°{M d ) and 



7 /n ■ 6 £X 

0(ar) = e d/2 f{ex)e 1 ^ 

be normalized in L 2 (R d ). Let a G C^(M d x IR d ) and T < oo. Then there is 
C < oo such that for all t G [0, T] 

i,af(t)0)- p ci (dxdp,t)a(x,p) 



< CVi, 



where p c \(dx dp,t) = (p c i o $ t )(dxdp) and p c \(dxdp) was defined in (pq 



(69) 



As in the preceding cases we can again translate Proposition [19] to the full 
dynamics using Theorem but we omit the corresponding statement this time. 
Proof. We apply Egorov and switch to standard quantization: 

(I a\m = (I (a o $*) W,£ ?) + 0(e) = ($, (a o $*)S,^) + 0(e) . 

We calculate for a G Cq° 

(a s '4)(|) = {2?r)- d e d/2 J dxdpa(ex,p)f{ex)e- ip < x ~ y/e) e lS(£x)/E 

= (2n)- d e- d/2 J dxdpa{x,p)f{x)e t{s{x) ' p - {x - y))/e 

= {2ir)- d e- d/2 J dzdpa{z + y,p)f{z + y)e l{S( - z+y) - p - z)/e . 

Stationary phase method (cf. Theorem 7.7.7 in fl with k = n + 1) yields that 

dzdpa(z + y,p)f(z + yy^^+yl-P-^^- 



(2ne) d a(y,VS(y))f(y)ei s ^ 



< Ce d+ i 



where the constant is uniform in y and depends on a via a sum of sup-norms of 
finite many partial derivatives. Hence 



((a o $*)S>4)(^) = e d,2 {a o VS(y))f(y)ei s ^ + e d/2 0(V^) 



, (a o = / dy ( a o $ ( )(y, V5(2/))/ 2 (y) + 0(y/i) 



and 



where we used that / G L 1 (H d ) 



□ 
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